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Abstract 


A model fot* the gLltltude dependence of the hot plesma perametere responsible 
for the electrostatic Oharging Of spOcect^aft has been developed. Based upon plasttia 
plasma orbit theory, the directed velocity is a function of the ambient magnetic 
field Huk density. A consequence of tnis approach is that \vhile the thermal velo- 
city distributions (assumed to be Maxwellian) of the plasma particles are independ- 
ent of the magnetic field strength (and hence altitude), the particle d: allies in- 
crease with ma^etic field strength. Thus, according to this model, While the 
equilibrium voltage is independent of altitude, the charging current density 
increases With decreasing altitude. Howevfer, the probability of such spacei^raft 
charging decreases with decreasing altitude. 


1. INTRODUCTION 

1 2 

Almost all of the published data and analyses of spacecraft charging * have 
been concerned with spacecraft In geosynchronous orbit (r « 6. 6 Rg). it is 
theoretically expected that the charging phenomenon can occur at other altitudes 

t 
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as W6U. J-or spacecraft In earth orbit at OthCr altitudes, the characteristics of 
I^h Charging becomes of practical concern, the purpose of this study ^as to 
develop an analytical model v^hich yielded the significant parameters Of the space- 
craft Charging phenomenon as functions of altitude above the earth. 

2. ALTltlJOE DEPENOESCE 

The major environment which has a string altitude dependence at geosjmchro- 
nous orbit is the earth-s magnetic field. If this field Is taken as approximately 
that of a dipole. US magnetic flUx density may be written 

B Jl + 3 sin^ X,. (1) 


where. 

B - the magnetic flux density at the surface of the earth at the magnetic 
° etiUator « 0.3 gauss = 3 X 10"® webers/m^: 

X = the magnetic latitude (degrees); 

r** = distance from the magnetic center of the earth: 

Rg = radius of the earth. 

It is necessary for r and Re to be in the same units, and that r > Re* 

It is vrell knovtn that a plasma can exist in a magnetic field only if the plasma 
energy density exceeds the magnetic field energy density. The energy ens y 

Uia medlutn aihtch conlalna the magnetic field. For free apace p - p 
H/m (MKS onlta). Thna for the dipole -approximated geomagnetic Held, 

(1*3 sin* i.,1 


^ j 1 + 3 sin 
3. 58 X 10"^ ^ j 


^M_1 job 
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Arty plasma which is able to exist In such a majgnatlc field must possess art enet*By 
density at least this large. 

The energy dehSlty of a two-component electrically neutral fully ionized rton- 
relativistic plasma may be writtert 


E. D.p = KE (directed) + KE (thermal) 


(3) 


where 


KE (directed) 
N 


1/2 (m^ + m_) V (dir)^ ; 


= 1/2 (th)2 + 1/2 m_ v_ (th)2. 

and m_ are the rest masses of the positively charged arid 
negatively charged plasma particles-, respectively, while v(dir) and v(th) are the 
directed and the thermal velocities of those particles. .. N is the spatial density of 
each type of particle. As long as the plasma moves as an entity, the + and - par- 
ticles will have the same directed velocities. In addition, if the pldsma is in 
thermal equilibrium, the -h and - particles will have the same average thermal 
leading to the relatibnship 


KE (thernial) 
N 

In these equations. 



This last relationship is only an approximation, as measurements in the hot plasmas 
responsible for spacecraft charging seem to show that the temperatures of the 

^ component (largely prOtohs) are about twice the temperatures of the - component 
(elebtrbns). 

The characteristics of the hot (spacecraft charging) plasmas at geosynchronous 
orbit show that the directed velocity generally lies between the thermal velocities 
of the electrons and the protons. Thus, the directed velocity is supersonic for the 
protons but subsonic for the electrons. SinCe m^ « ld36 m , the energy density 
for such plasma is essentially all due to the directed motlon’of the protons. To 
illustrate this, for a plasma with a 5 keV temperature v^(th) a 10® cm/sec while 
V.(th) a 4.3 X 10 cm/sec. Taking v(dlr) as 6 X 10® cm/sec (an intermediate 
value) yields directed energies of E^ a isO keV and E. a 0. 1 keV. Only if v(dlr) 
is not greater than V^(th) is the plasma energy density not largely due to the 
directed motion of the heavier (+) component. 

With E. D. p a 1/2 m^ v(dtr)2, it is possible to estimate how low in altitude a 
hot Plasma moving radially toward the earth can get before Its energy density Is 
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not sufficient to prevont the geomagnetiO field from tearing It apart. Equating 
D. g and E. D. yields 


!•» Rj; 



1 + 3 Sln2 

l^vntdTrTj 



4S) 


The results of this calculation are shown graphically in Figure 1 for ■ 0*^ 
(magnetic equator). These results are an approximation because of many factors, 
but due to the steep radial gradient of the energy density of the geomagnetic field. 



2 do 



the approxirtations pi-oduee ohljr relatively small perturbations In the results. It 
Is seen that a moderately dense (> lO partleles/em^) plasma with a reasonably 
large directed velocity (^3 X 16® cm/sec) can reach the 12 hr circular orbit 
(r e 4. 15 11^), an Increase In either particle density or directed velocity 

can result in spacecraft charging at even lower altitude earth orbits. 


3. modeling 

There are various approximate ways of modeling a plasma to make plasma 
problems mathematically tractable. Each model deals with that portion of plasma 
properties which are most pertinent to the problem at hand while suppressing or 
ignoring less relevant properties. Plasma orbit theory is among those approximate 
models often used in situations in which the motions of individual plasma particles 
are important. (Hydrodynamic theory is often used for situations in which large- 
scale plasma oscillations are more important than individual particle motion. ) 

The basis of orbit theory is the Conservation of the angular momentum of individual 
particles about an axis (for example, the direction of the magnetic flux lines). 
Motion of the piasma particles parallel to the lines of magnetic flux is not affected 
(to a first approximation). Thus, it is possible to separate the thermal and the 
directed motions of the plasma particles by considering them to be parallel and 
perpendicular to the magnetic field, respectively. This simplified model results 
in a decrease in the directed velocity as the distance from the earth decreases 
(for example, as the magnetic field strength increases) while leaving the thermal 
velocity unchanged, it is obvious that this model neglects the compressional 
heating of the plasma, which certainly is a major factor in the production of the hot 
plasma in the first place. However, attempts to incorporate compressional heat- 
ing resulted in a model in which the plasma energy density increases as rapidly 
as the geomagnetic field energy density, with the consequence that a plasma able to 
exist at any altitude could theoretically reach any other (lower) altitude. This is 
clearly at variance With observational data (spacecraft charging anomalies Have 
not been reported for low -altitude earth orbits). 

Based upon the simplified application of plasma orbit theory. It is possible 

to write 

v(dlr) = cos a 

where a the pitch angle (angle bet^veen the velocity vector* and the dlrectlori of 
the geomagnetic field), this Is similar to the ecjuatlort of motion for Van Allen belt 
particles, the second Invariant for such particles leads to the relationship 
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By eomblAtng thbse last two equations, onb can obtain 


'^dlr = '^th 


V 


fi 

B • 

mak 


(8) 


According to this model the average particle thermal Velocity (plSsma tempera- 
ture) does not change Vvith altitude while the directed velocity gradually decreases, 
reaching sero when B = point, the plasma energy density has been 

reduced to the energy density Of the geomagnetic field. 

The negative voltage to which a body immersed in a plasma of a given tempera- 
ture will change is 76 times the average electron kinetic energy (e!^) in electron 
Volts. The negative voltage is due to the fact that the electrons have much higher 
thermal velocities than the protons, and therefore impact the spacecraft surface 
much more often. 

The factor of 3. 76 is due to the fact that the negative current to an unillumi- 
nated spacecraft Surface varies exponentially with spacecraft Voltage (to a first 
approximation) While the positive current varies linearly With voltage. This is a 
consequence of the fact that the electron motion is subsonic while the proton motion 
is. S upe r s onic . Mat he ma tic ally. 



where and are the surface current densities as a function of Spacecraft 
voltage (V), and and are those current densities (including the effects of 
secondaries) when V = O. V^p and are the average proton and electron kinetic 
energies In electron volts, respectively. Since m_^^ « 1836 m^, ^ 18^6 

« 42. 8 J^p. The value of V will increase until = J^, assuming no electrical 
discharges take place. Equating J ^ to and solving for V yields 





( 11 ) 


Silice th6 secondary currSAt components Vary with spacecraft voltage differently 
than the primary currents do, this factor of j.76 cati be appreciably different in 
many sttuatlohs. HOVrever, sihce the particles responsible for the charging have 
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art rtppt'dximateiy Maxwelllart energy distribution, It Is understandable that a plasma 
with an average cleetron temperati-re of V^g (volts) woiild charge a spacecraft to 
voltages V > Vgg. 

4. CUARGINU-PARAMBTI^OS 

Since, according to this model, the particle thermal Velocities are unchanged 
by the directed motion of the plasma toward the earth, the equilibrium Voltage 
which a spacecraft surface Will reach in the plasma will not depend upon altitude. 

In the absence of electrical discharges, a spacecraft in a 12 hr circular orbit 
(r ^ 4. 15 Rj.) will ’herefore theoretically charge to the same potential as one in 
a 24 hr geosynchronous Orbit <r s 6. 6 Rg). However, the surface current densi- 
ties (which determine how fast the spacecraft will charge) will be a function of 
altitude. The initial primary surface current density as a function of altitude may 
be estimated by calculating the limiting plasma particle density of the plasma jUst 
as its directed motion has ceased. Since the average particle energy (and the 
particle energy distribution) remains Unchanged according to this model, the par- 
ticle density must increase linearly with the energy density of the geomagnetic field. 

■the Calculated numbers derived based Upon this model are shown in Table 1. 

At each altitude r the average energy density of the gfeomaghCtic field (dipole 
approximation) was taken as the average thermal energy density of a stationary 
plasma at that altitude. Assuming half of this thermal energy density was dUe to 
the electrons, the product • N_ was obtained. For each value of limiting space- 
craft potential V, the average electron energy was obtained by dividing by 3. 76. 

The electron density (Cm‘3) was the quotient of E • divided by (E sE here). 
The initial primary electron rrent density (J_) was obtained from the equation 

J_ = N_ * q • V . (12) 

The results listed in Table 1 are also shown granhlcally in Figure 2. It is 
seen that the calaulated initial current density increases rapidly as orbit altitude 
decreases. This altitude dependence decreases the importance of the photoelectric 
current density (Which has an altitude -independent value < i nA/cm^) at lower 
altitudes. Another consequence of this altitude dependence is that if the spacecraft 
surface cannot withstand the equilibrium potential to Which the plasma can charge 
it, the rate at which electrical breakdowns (discharges) occur will be much greater 
at lower altitudes. 

While the calculated surface current densities are larger at low altitudes than 
at geosynchronous altitude, the probability that a spacecraft will encounter the hot 
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1 

V. ./cm*) 

V • 3 IteV 
ie- 1.55 ItaV 
■ - a.s^io* 

V - lO keV 
Ee- 2.88 kbV 
0 - 8.08xlO< 

V - Is kbV 
Ee- 5.99 keV 
9 • 7.45x10* 

V - 20 koV 
Ee- 5.32 keV 

V - 8.6x10* 

V - 25 keV 
Ee- 6.65 keV 
9 - 9.62*10* 

6.6 

6.5 

k. - 4.34 

2.44 

1.63 

1.22 

0.97 

(24 hr) 


J- - 3.36 

2.3d 

1.94 

1.68 

1.50 

6 

12.5 

M_ - 9.40 

4.70 

3.13 

2.35 

1.88 



J. ■ 6.47 

4.57 

3.74 

3.23 

2.89 

5.5 

20 

N. • 15.0 

7.52 

5.00 

3.76 

3.01 



d. - 10.35 

7.3d 

5.96 

5.17 

4.63 

5 

36 

H_ - 22.1 

13.5 

9.00 

6.80 

5.40 



J. * 18.6 

13.2 

10.7 

9.29 

8.31 

4.5 

67,5 

K. - 50.8 

25.4 

16.9 

12.7 

10.2 



J_ - 34.9 

24.7 

20.1 

17.5 

15.6 

4.15 

lid 

a. - 82.7 

41.4 

27.6 

20.7 

16,6 

{\2 hr) 


j. - S8.8 

40.2 

32.8 

28.5 

25.5 

3.5 

313 

H_ - 235.0 

117.5 

78.3 

58.8 

47.0 



J_ - 161.5 

114.3 

93.3 

80.8 

72.3 


V In (im/a6fei A. ib p&ttlbles/cm’; J_ lb namps/cm^. 


Table 1. Calculated Partible Densities iuid initial Pfitnary Current Densities as 
Functions of Altitude = 0®) 

plp-dma at these Ibwer altitudes is appreciably less. It would be theoretically 
possible to calculate an altitude -dependent probability if data arere available con- 
cerning the pt^jbability of encountering a given plasma energy density at geosyn- 
chronous altitude. Such data are becoming available. However, the measured 
L -dependence of the probability of encoiiAtering >15 volts on an antenna has been 
measured by lMP-6^ (see Figure 3). This spacecraft was launched 31 March I9tl 
into a 26.7° elliptical orbit with an apogee of 32.4 Rg and a perigee of 243 km. 

As the data shows, the probability of encountering hot (Or at least warm) plasma 
at 4. 15 Rg is approximately an order of magnitude less than at 6. 6 Rg. While the 
inclination of the IMP-6 orbit makes analysis difficult, the data suggest that the 
spacecraft charging phehomenum may be more comtnon near geosynchronous 
altitude than at any other. 
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Figure Prubabtlity of Measuring > 15 Volts on IMP>6 Antenna 
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